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We construct a family of rotationally invariant, local, S=l/2 Klein Hamiltonians on various lattices 
that exhibit ground state manifolds spanned by nearest-neighbor valence bond states. We show 
that with selected perturbations such models can be driven into phases modeled by well understood 
quantum dimer models on the corresponding lattices. Specifically, we show that the perturbation 
procedure is arbitrarily well controlled by a new parameter which is the extent of decoration of 
the reference lattice. This strategy leads to Hamiltonians that exhibit i) Z2 RVB phases in two 
dimensions, ii) U(l) RVB phases with a gapless "photon" in three dimensions, and iii) a Cantor 
deconfined region in two dimensions. We also construct two models on the pyrochlore lattice, one 
model exhibiting a Z2 RVB phase and the other a U(l) RVB phase. 



I. INTRODUCTION 

Just over thirty years ago Anderson^ introduced the 
resonating valence bond (RVB) state as an alternative to 
Neel ordering in antiferromagnets with strong quantum 
fluctuations. In essence, he proposed that on a suffi- 
ciently frustrated lattice an S = 1/2 system would ex- 
hibit a disordered state at T = 0, which would be cap- 
tured by a wavefunction of the form 



= 5> 



(1.1) 



where |c) is a configuration of singlet pairings of spins 
or valence bonds (Fig . For sufficiently short-ranged 
valence bonds this describes, in contrast to the Neel state, 
a state with short-ranged spin correlations. 
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FIG. 1: Sample valence bond configurations. The thick lines 
represent singlet pairings. The original formulation allowed 
for bonds of arbitrary length (left). Considerable progress has 
been made by restricting to configurations with only nearest- 
neighbor valence bonds (right). 



The discovery of the cuprates and the suggestion that 
their superconductivity could be traced to RVB physics 2 
greatly energized the elucidation of the RVB idea and 
by now a rather complete understanding of its inter- 
nal logic has emerged. In modern parlance, an RVB 



phase is a topological phase, characterized by excita- 
tions with fractional quantum numbers and a low energy 
gauge structure which mediates topological interactions 
among the excitations^. The excitations include spinous, 
the S = 1/2 excitations produced by breaking a valence 
bond, as well as collective excitations within the valence 
bond manifold (see figure[5J). The sr-RVB 4 (short ranged 
RVB) with short-ranged bonds and gapped spinons will 
be our concern in this paper. A post-cuprate version with 
longer-ranged bonds and gapless spinons^ has also been 
the subject of recent progress £ Readers familiar with one 
dimensional lore will note that the short-ranged and long- 
ranged RVBs generalize the physics, respectively, of the 
Majumdar-Ghosh chain 7 and the Bethe chain 8 to higher 
dimensions. 



FIG. 2: Typical excitations of an RVB liquid. The left fig- 
ure depicts a pair of spinons, fractionalized excitations with 
S = 1/2, formed by breaking a valence bond. Spinons are 
unpaired spins which move in the RVB liquid background; 
the interaction between spinons depends on the lattice geom- 
etry. The right figure depicts a vison, which is an excitation 
within the valence bond subspace. If we consider an RVB 
liquid which is an equal amplitude superposition of nearest- 
neighbor valence bond states, then the vison is defined by the 
wavefunction \ip v ) = X* ( — -0 l c ) where N c is the number of 
bonds which cross the dashed line shown in the figure. Clearly 
this state is orthogonal to the RVB state, \tp) = ^ c \c). The 
interaction between spinons and visons is discussed in detail 
in Ref. [| 
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The important progress that we have described has 
been kinematical. It has not directly answered the ques- 
tion of realizing RVB phases for actual Hamiltonians. 
The original proposal was made for the nearest-neighbor 
Heisenberg model on the triangular lattice but that is 
now generally believed to exhibit weak Neel ordeA To 
make progress on the dynamical front, Rokhsar and 
Kivelsonifi introduced the quantum dimer model (QDM) 
which assumes that the low energy dynamics is dom- 
inated by valence bond configurations of short range 
which are taken to be nearest-neighbor in the versions 
studied to date. Such configurations are labelled by 
dimer coverings of the lattice at issue and the quantum 
dimer Hamiltonian acts in a Hilbert space spanned by 
such coverings. The program of studying the simplest 
dimer models has been rather fruitful. It is now clear 
that Zi RVB phases may arise on non-bipartite lattices 
in d > while bipartite lattices in d > 2 give rise to 

U(l) RVB phases that exhibit a gapless "photon"^ii4. In 
addition a variety of crystalline phases have been identi- 
fied, most notably a Cantor deconfined region^ of inter- 
leaved commensurate and incommensurate valence bond 
crystals on bipartite lattices in d = 2. 

The next order of business then, is to find rotation- 
ally invariant, local, spin Hamiltonians that are accu- 
rately described by these well understood dimer mod- 
els. This is the problem that we solve in this paper thus 
completing a program initiated by Chayes, Chayes and 
Kivelsoni^. The strategy we follow is that of constructing 
Klein Hamiltoniansii with large energy scales which se- 
lect nearest neighbor valence bond states as their ground 
states, separated by a gap from excited states. We then 
lift this degeneracy by the inclusion of perturbations that 
precisely mimic the terms in the quantum dimer models 
of interest. To control this procedure, with its difficul- 
ties stemming from the non-orthogonality of the valence 
bond basis, we introduce a parameter which is the extent 
of decoration of the reference lattice. By suitably tun- 
ing this parameter we are able to make our dimer model 
realizations arbitrarily accurate. An elegant feature of 
this limit is that it enables us to establish the existence 
of a gap about the nearest-neighbor valence bond mani- 
fold and to discuss the states above this gap in terms of 
"microscopic" spinons whose meaning will become clear 
below. 

We note that our tuning procedure occurs within the 
space of SU(2) invariant Hamiltonians. This is in con- 
trast to approaches involving enlarging the symmetry 
group to Sp(iV) or SU(iV) and studying the large N 
limi1*i&; in these cases, the applicability of results to Sp(l) 
= SU(2) is not obvious. 

This is also a good place to note that there is a con- 
siderable body of work on variational 19 and finite-size 
studies of some of the two-dimensional phases discussed 
in this paper, e.g. the early finite-size study of a multiple- 
spin Hamiltonian on the triangular lattice that adduced 
evidence for a topologically ordered phased A compre- 
hensive review of such work is given in Ref. EH This 



work is complementary to ours as it deals with some- 
what simpler Hamiltonians but is unable to access the 
thermodynamic limit in a controlled fashion. We also 
note that, there is a large and growing literature on more 
general models with topological phases which we skip in 
our focus on S = 1/2 spin systems; this work also finds 
inspiration from the proposal that a quantum computer 
may be robustly created from a topological phased Fi- 
nally, we note that there are encouraging reports of spin 
liquids in experimental systems^^ 

In the rest of the paper we give details of our construc- 
tions. We begin with a quick review of quantum dimer 
models and the known results on their phase diagrams in 
Section II. In Section III we explain our strategy with the 
honeycomb lattice serving as an example; this realizes the 
physics of bipartite dimer models in d = 2. In Section IV 
we show how the physics of non-bipartite dimer models 
in d = 2 and bipartite and non-bipartite dimer models 
in d > 2 may be obtained from spin models. In Section 
V we discuss two spin models on the non-bipartite py- 
rochlore lattice, one exhibiting a Z2 RVB phase and the 
other a U(l) RVB phase. We conclude with a summary 
(Section VI) and a set of appendices that contain some 
technical material. 

II. QUANTUM DIMER MODELS 

For a system of spins on a lattice A, a (nearest- 
neighbor) valence bond state is a product wavefunction 
of the form * = ipij where ipy = ^W'JV'j - V'j'V'j) 
and the product is over nearest-neighbor pairs (i, j). The 
product is defined so each spin forms a singlet with ex- 
actly one of its neighbors. Each valence bond state corre- 
sponds to a hard-core dimer covering of A where a dimer 
connecting two sites corresponds to a singlet bond be- 
tween the respective spins. Valence bond states are not 
orthogonal but the overlap between two arbitrary states 
is exponentially small in the length of closed loops ob- 
tained by superposing them. This suggests that they are 
linearly independent on sufficiently open lattices and in- 
deed there are proofs for manyi£ and numerical evidence 
that this is so even on the triangular lattice^. 

The identification with dimer coverings suggests that 
any low energy dynamics restricted to the valence bond 
manifold can be represented by a quantum dimer Hamil- 
tonian acting on orthogonal dimer states. The simplest 
such Hamiltonian on the square lattice, written down by 
Rokhsar and KivelsoniS, has the pictorial form 

H qdm = 2-t(|l I )(n I +h.c.)+v(| r I ><I 1 1 +|~><~ I ) 

(2.1) 

where t and v are positive constants and the sum is over 
all possible square plaquettes. Evidently, this can be sup- 
plemented by kinetic energy terms which act on longer 
loops and potential energy terms which count more com- 
plicated dimer motifs^ 
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The passage from valence bonds to dimers, however, 
has to contend with two complications. One is that one 
needs to choose a phase convention for the valence bonds, 
which is subject to restrictions on what signs one can 
obtain for various couplings in the dimcr Hilbert space. 
The other, already alluded to, is the lack of orthogonality 
of the valence bond states which makes the transcription 
from a spin model to a dimer model non-trivial. We 
will deal with both problems later in the paper; here we 
merely wish to alert the reader to their existence. 

An important property of the quantum dimer Hamil- 
tonian is the existence of the "Rokhsar-Kivelson point" 
(RK point) t = v, where any equal amplitude superpo- 
sition of all dimer coverings connected by the operation 
of the kinetic energy is a ground state. To see this, note 
that for every flippable plaquette, the second term gives 
a penalty v while the first term gives at most a benefit 
of —t. Nonflippable plaquettes are destroyed by Hqdm- 
This gives a lower bound for the ground state energy: 
Eq > min{0, Np(v — t)}, where Np is the number of 
plaquettes in the lattice. The equal amplitude state has 
energy (nfi)(v — t), where (rifi) is the average number of 
flippable plaquettes in the state. At v = i, this saturates 
the lower bound and, since the equal amplitude state is 
an eigenstate of Hqdm ( a t v — t), we may conclude that 
it is the ground state when v — t. 

Thus the ground state correlations at the RK point 
reduce to those of solvable classical dimer models. Ad- 
ditionally, the infinite temperature static correlations of 
QDMs also reduce to those of the same classical models. 
These features, along with the additional one that Hamil- 
tonians of the form i|2.1[l can be simulated by Monte- 
Carlo without any sign problems, have been crucial to 
making progress in determining the phase diagrams of 
the quantum models. 

As a consequence of this progress we now know that: 

i) QDMs on bipartite lattices in d = 2 do not exhibit 
an RVB phase. The equal amplitude state present 
at the RK point, v/t = 1, does not extend into a 
phase. As v/t increases, the system generically passes 
through a sequence of interleaved commensurate and 
incommensurate crystalline phases before reaching the 
staggered valence bond solid (VBS) phase. These inter- 
mediate phases, whose measure approaches unity near 
the RK point, turn out to have deconfined monomers, 
a phenomenon coined Cantor deconfinemenljA^. As 
v/t decreases from unity, the system passes through 
a plaquette phased before undergoing a first-order 
transition to a columnar VBS. 

ii) QDMs on non-bipartite lattices in d = 2 may exhibit 
RVB phases. These are Z 2 RVB phases captured by 
a purely topological BF theory^. For the triangular 
lattice, it has been shownAi that for v/t > 1, the system 
is in a staggered VBS; for v/t < 1, there is a deconfined 
RVB liquid phase. As v/t is further reduced, there are 
probably a small number of VBS phases culminating in 
the columnar state. 

iii) QDMs on non-bipartite lattices in d = 3 and higher 



also exhibit a Z 2 RVB phased. 

iv) QDMs on bipartite lattices in d = 3 and higher 
exhibit a U(l) RVB phase with a gapless, linearly 
dispersing transverse mode, the "photon"—. 

We now turn to the task of constructing spin models 
whose low energy dynamics is precisely captured by these 
dimer models. We begin, for pedagogical specificity, with 
the honeycomb lattice. 



III. HONEYCOMB LATTICE: BIPARTITE 
PHYSICS IN d — 2 

Our strategy for realizing dimer models proceeds in 
three steps. First, we construct, following Klein, a lo- 
cal spin Hamiltonian that has valence bond states as its 
ground states. Next we perturb it to obtain a QDM. Fi- 
nally, we decorate the lattice to simplify the QDM to the 
well studied form 1)2.1(1 . In the high decoration limit, we 
show the existence of a gap and give a description of the 
spectrum in terms of spinons. 



A. Klein model 

The basic idea of the Klein model consists of consider- 
ing a cluster of z sites (typically, a spin and its z— 1 neigh- 
bors) and deterring, via an energy penalty, this cluster 
from having maximal total spin Stot = z/2. If two of the 
spins in the cluster form a singlet bond, this condition is 
satisfied. This is why Klein Hamiltonians naturally lead 
to valence bond ground states. 

In particular, for a system of spins si on a lattice A, a 
Klein Hamiltonian is a sum of projection operators 
defined as follows. For each site i, consider the neighbor- 
hood of spins N(i) consisting of the spin at site i and 
its (zi — 1) nearest neighbors. Let fjv(i) be the projector 
onto the highest total spin state of the cluster. The Klein 
Hamiltonian is formally given by the expression: 



with total spin of cluster Af(i) given by: 



(3.1) 



(3.2) 



We may write fV(i) in terms of this operator. For exam- 
ple, if Zi is even, then: 



n [* 

L=0 



'AT(i) 



L(L + 1] 



(3.3) 



The total spin of this cluster will take values from 0, 
1,. . . ,(zi/2) — l,Zi/2. The factors of this product are op- 
erators which sequentially annihilate all but the highest 
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spin sector. The form of the operator implies that larger 
clusters involve higher-order spin interactions- but they 
always remain local. 

If the constants Ci in Eq. I3.3l are chosen to be positive, 
then Hk will have non- negative eigenvalues. By con- 
struction, valence bond coverings are zero-energy ground 
states of Hk- 

For the honeycomb lattice, Chayes et ali& have already 
written down the explicit form of the Klein Hamiltonian 
in terms of spin operators. Their expression, up to unim- 
portant overall constants, is: 

H = ^ Si- Sj + - g i ■ 

2 ' 

+ 5 zJC^'^OO 5 * - g r) (3.4) 

ijkr 

The first and second terms are over nearest and next 
nearest neighbors respectively. The third term is over 
quartets ijkr where i and j are nearest neighbors; k is 
a neighbor of i different from j ; and r is a neighbor of j 
different from i. A striking feature of this Hamiltonian 
is that the leading term is the usual Heisenberg antifer- 
romagnet. 

B. Perturbations 

We will now perturb the Klein Hamiltonian to obtain a 
QDM with dynamics. In doing so we will use the overlap 
expansion invented by Rokhsar and KivelsoniS. . which 
is predicated on the linear independence of the valence 
bond states. That valence bond states on the honeycomb 
lattice are linearly independent was proved in Ref. IT^ 

For the purpose of obtaining the dimer kinetic energy, 
it is sufficient to consider including just an additional 
nearest neighbor interaction^, 

SH = ^ s, ■ s 3 (3.5) 

To first order in degenerate perturbation theory, we may 
write this as an effective operator on the valence bond 
states. First, we define an orthonormal basis set {\a)} 
in terms of the linearly independent valence bond states 

my- 

R = E(^ 1/2 )^K> (3-6) 

i 

Here Sij = (i\j) is the overlap matrix element between 
valence bond states \i) and \j). The magnitude of the 
overlap of two valence bond states may be determined by 
overlaying the two configurations forming what is called 
the transition grapk 1 ^. The construction is described in 
Fig. El As shown in Fig. |3 the transition graph consists 
of double bonds, where the two states have a bond in 
common, and closed loops of varying (even) lengths. The 



magnitude of the overlap S!y is given by 2 Nl J]; x L ' where 
Ni is the number of loops; the product is over all such 
loops; Li is the length of loop i; and x — The sign of 
Sij depends on how we choose to orient the bonds. By 
orientation, we refer to the fact that a bond between sites 
1 and 2 may be interpreted as the singlet bond tpis = 
-j^{i>\i>2 — V'iV'I) or as V"2i = — V'12- The key idea of 
the overlap expansion is to treat X (IS cl small expansion 
parameter. 




FIG. 3: Transition graph construction for two valence bond 
coverings of the square lattice (the construction for other 
lattices is similar). The dots are the lattice sites and the 
thick lines denote that the two sites form a singlet bond. 
The singlet orientation may be specified, for example, by 
having the bonds point from the gray sites to the white 
sites. The transition graph is formed by overlaying the two 
configurations resulting in a graph (lower) containing dou- 
ble bonds and closed loops of varying (even) lengths. In 
the above example, there are two double bonds and three 
loops of lengths 4, 6, and 16. The magnitude of the over- 
lap between the two valence bond coverings is then given by 
\S\ = 2 3 (l/ % /2) 4+6+16 = (i) 10 = T -L. Thus, while the over- 
lap between two arbitrary valence bond coverings is never 
zero, it is usually a small number. This is the basis for the 
overlap expansion discussed in the text. 



We may specify the matrix elements of our effective 
operator in terms of the {\a}} basis: 

H a0 = {S-^SHS- 1 ' 2 )^ (3.7) 
= Y,{S-V*U(i\8H\f){S-V*) sf > (3.8) 

ij 

If either state \i) or \j) contains the bond (12), then (i\si ■ 
s*2 |j) = — 11 neither \i) nor \j) contains the bond 

(12), then we have a non-zero matrix element only if spins 
1 and 2 are members of the same loop in the transition 
graph. If that is the case, then one may show that (i\si ■ 
82\j) = (^l)"^ 2 (T|)(«|j) where n is the length of the 
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loop for the case where spins 1 and 2 are separated by 
an even (odd) number of sites. 

We now specialize to the honeycomb lattice. As 
sketched in the Appendix, we may orient the bonds on 
the honeycomb lattice so that for any two states differing 
by a (minimal) length 6 loop, the overlap has a pos- 
itive sign. On the honeycomb lattice, the sign is a matter 
of convention and we could have chosen the negative sign. 
Given the positive sign convention, we conclude that our 
matrix element is given by: 

(i\SH\j) = -^5 l3 - \{2x G )o %3 + 0(x 10 ) (3.9) 

where is the number of bonds (half the number of 
sites) and Ojj is a matrix that is 1 if states \i) and \j) 
differ by a length 6 loop and zero otherwise. We may 
also expand the overlap matrix: 

Sij = 6 tj +2.x 6 O y - +0{x 10 ) (3.10) 
(S~ 1/; %- = % - x 6 O l3 + 0(x 8 ) (3.11) 

Comparing the previous line with Eq. 13.61 we see that 
within the overlap expansion, each |a) has a largest com- 
ponent corresponding to a unique valence bond state. 
Therefore, we will refer to the orthogonal set {\a}} as the 
set of dimer states corresponding to the valence bond cov- 
erings. In writing our effective operator, we absorb the 
leading term involving the number of dimers rid times the 
unit operator into our definition of zero energy. 

H af3 = ^[^-a; 6 O m + 0(x 8 )] 

ij 

x [-^(2x 6 )o l3 + 0(x 10 )} 

x [S j0 - x 6 O jP + 0(x s )} (3.12) 

« ~(2x 6 )o a0 + O(x 10 ) (3.13) 

We conclude that the leading term in the overlap expan- 
sion is an operator with nonzero matrix elements only 
between dimer states differing by a minimal length 6 
loop. All of these nonzero elements have the same value 
—t = — (3/4)(2a; 6 ). Thus we have obtained the kinetic en- 
ergy operator in the quantum dimer model (QDM). Note 
that we can conclude this only because we were able to 
define a bond orientation convention such that all mini- 
mal overlaps come with the same sign. Otherwise some 
off-diagonal terms would have energy t and the energy 
bound arguments which we gave previously to conclude 
that there is an RVB state at the RK point would no 
longer hold. 

To obtain the potential energy term as the leading or- 
der effect, we need a more complicated interaction which 
we take to be: 

SH = J~^^ Si ■ Sj + 

(ij) 



v'Ylifa • s 2 )(s 3 • s 4 )(s 5 • s 6 ) 
o 

+ (s2 • s 3 )(s 4 • s 5 )(s 6 • si)^) (3-14) 

where the first sum is over nearest neighbors and the 
second sum is over elementary plaquettes (see Fig. 0J. 




FIG. 4: The elementary plaquette of the honeycomb lattice. 

The first term of Eq. I3~T1 gives the QDM kinetic 
energy. Similarly, consider the operator S12S34S56 = 
(si • S2XS3 • S4XS5 • sq). If valence bond states \i) 
and |j) contain the bonds (12), (34), and (56), then 
(*|si2S34S56|j) = (—j) 3 (i\j)- If state \i) contains all three 
bonds, then the diagonal matrix element (?|si2S34S56K) 
is (— |) 3 . If state \i) is missing one or more bonds, then 
the diagonal matrix element is zero unless |z) contains 
the three complementary bonds (23), (45), and (61). In 
this case, the expectation of the operator equals (— j ) 3 x 6 , 
which is higher order in the overlap expansion. It may 
be shown that off-diagonal matrix elements evaluate to 
a term proportional to the overlap of the states, the pro- 
portionality constant being of order unity. These results 
imply that the matrix element of the ring interaction be- 
tween two valence bond states is given by: 

(-^fvn fUl 5 l3 +0{vx 6 ) (3.15) 

where rifij is the number of flippable hexagonal plaque- 
ttes in configuration |z). We may write this in terms of 
dimer states, as discussed above. Absorbing numerical 
factors into the constants J and v, we arrive at our ef- 
fective dimer Hamiltonian: 

H a p — — Jx 6 o a p + vrifi ia S a /3 + 0(vx e + Jx w ) 

= —tO a p + vnfi ia 5 a fi + 0(vx e +tx i ) (3.16) 

where t = Jx e and n/i, Q is the number of flippable pla- 
quettes in the valence bond state corresponding to dimer 
state a. If t and v are of order unity, then the higher or- 
der terms will be small compared to the first two terms, 
which act on our dimer states (which are really spin 
states) in a manner analogous to the QDM kinetic and 
potential energy operators on usual dimer states. 
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For the actual problem at hand, x — (1/V2) is less 
than 1 but is by no means tiny. Hence the neglect of 
other terms induced by our perturbations is not obviously 
justified. While we do not need them to be zero, we 
do need them to be weak enough perturbations so the 
analysis of Ref. is justified. 

What we do learn from the overlap expansion result 
(|3.1(il) is that the non-orthogonality is a much smaller 
problem on more open lattices which involve large loops. 
While the honeycomb is a good candidate on this score, 
to put the issue beyond doubt we now consider a deco- 
rated version of the lattice. 



C. Decoration scheme 

In this section, we propose a modification to our earlier 
arguments which makes the overlap expansion essentially 
exact. Consider the decorated honeycomb lattice shown 
in Fig. [5] where we insert N (an even integer) sites be- 
tween neighboring sites of the usual honeycomb lattice. 
The dimer structure of this lattice, including the number 
of dimer states, is exactly the same as before except that 
having a dimer between sites 1 and 2 corresponds to a 
chain of (N + 2)/2 dimers beginning at site 1 and ending 
at site 2. Not having a dimer between sites 1 and 2 cor- 
responds to having a chain of N /2 dimers beginning at 
site b\ and ending at site a^. The Klein Hamiltonian is 
correspondingly modified by including Klein projectors 
for the added sites. 

Majumdar and Ghosbi showed that the valence bond 
state is the only ground state of the Klein Hamiltonian 
for a one-dimensional spin chain with an even number 
of spins. Therefore, the conclusions regarding the Klein 
model on the honeycomb lattice (linear independence of 
valence bond states, valence bond states span the ground 
state manifold, etc) carry over directly to the decorated 
honeycomb lattice. 

While decorating does not introduce any new technical 
problems, there is a significant technical advantage with 
respect to the overlap expansion. The smallest two loops 
on the hexagonal lattice are length 6 and 10 from which 
we obtained that the relative orders of the leading and 
error terms in the overlap expansion were x 6 and a; 10 . 
The smallest loops on the decorated hexagonal lattice 
have lengths 6(N + 1) and 10 (iV + 1). Repeating the 
previous analysis, we will find that the leading and error 
terms in the overlap expansion are x 6 ( N+1 ^ and x 10 ^ N+1 \ 
The ratio of error term to leading term has improved 
from x 4 to x 4N . In the large N limit, the error term is 
"rigorously" negligible but we propose that even fairly 
small values of N may suffice to capture the qualitative 
features of the large N limit. 

While we have added complexity to the lattice, we do 
not have to increase the order of the spin interaction. 
Consider the following as a perturbation to the decorated 




FIG. 5: The decorated honeycomb lattice where N (an even 
integer) two-fold sites are inserted between the old sites. This 
drawing shows N = 4. The labels ai and b\ designate the 
first counterclockwise and first clockwise neighbor of spin 1 
where clockwise is with respect to the loop 123456. For the 
undecorated case, ai and b\ are just sites 6 and 2. 



honeycomb lattice Klein model: 
SH = J Si ■ Sj + 

o 

+ {si ■ s ai )(s 3 ■ s a3 )(s 5 ■ s 05 )J (3.17) 

The first term is a nearest neighbor interaction over all 
spins while the second term is over all elementary pla- 
quettes, such as the one in Fig. [SJ A 6 spin interaction is 
sufficient, even though we have many more spins in the 
loop, because having a (l&i) bond automatically implies 
the other bonds in the chain connecting 1 and 2. Our 
previous analysis carries over to the present case and we 
conclude: 

H aj3 = -Jx 6(N+1) O a0 + vn f i, a 5 a f3 
+ 0(vx^ N+1 ^+ Jx w ( N+ V) 
= -tO aP + vn f i, a S a0 + 0(vx^ N+1 ^ + tx AN ) 

(3.18) 

where t — Jx^ N+1 ^ and otherwise the notation is the 
same. 

Clearly, by decorating enough we can make the matrix 
elements beyond the dimer model arbitarily small and 
thus realize the physics, including Cantor deconfincmcnt, 
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present in generic, weak perturbations of the honeycomb 
QDM. 



1. Spinous 

In the highly decorated limit, one may show that 
nearest-neighbor valence bond states are separated by a 
finite gap from the excited states of the Klein model. In 
this limit, we are connecting a set of Majumdar- Ghosh 7 
(MG) chains into a two-dimensional network. We may 
describe the excited states of our system in terms of the 
well studied spinon defects of the MG chains, which are 
widely believed to be gappedS^ii^i. Here we give an 
outline of our argument and relegate technical details to 
Appendix [5] 

In Appendix [5] we consider what happens when we 
put these chains together for different values of a tunable 
parameter in our model: the ratio of Klein scales (the 
coefficient C, in Eq. I3.3|l of the Klein projectors of the 
decorated and reference sites. For large values of this ra- 
tio, the excited states are represented by "microscopic" 
spinons localized on the reference sites. For small val- 
ues of the Klein ratio, the excited states are extended 
and may be interpreted as MG spinons scattering at the 
vertices. There is a first order transition between these 
limits. In both limits, there is a gap between the VB 
manifold and the spinon states, as depicted in Fig. |SJ 

Spinons are the natural excited states (outside of the 
VB manifold) in the high decoration limit. For an un- 
perturbed Klein model, the VB manifold is degenerate so 
these excitations are mobile. The next question is what 
happens when the degeneracy of the ground state man- 
ifold is lifted. We argue that this has a small but vital 
effect on the spinon dynamics. At the RK point and in 
liquid phases, we expect the spinons to be deconfincd. 
In the crystalline phases, we may consider a pair of test 
spinons, holding one member fixed and considering the 
quantum mechanics of the other. If the wavefunction 
of the non- fixed spinon has spatial extent L, this would 
have an energy cost of order e c L 2 where e c is the en- 
ergy cost per unit area of scrambling the crystalline back- 
ground. The L-dependence of the kinetic energy varies as 
1/L 2 . The implication is that while e c is a much smaller 
scale than the spinon gap, spinons moving in a crystalline 
background are confined at sufficiently long length scales. 



D. Square lattice 

The square lattice is another two-dimensional bipartite 
lattice for which (nearest-neighbor) valence bond states 
are linearly independent^. We may orient the bonds on 
the square lattice so that two states differing by a (min- 
imal) length 4 loop, have positive overlap. The problem 
with applying our approach to the bare square lattice is 
that the Klein model has ground states outside of the va- 
lence bond manifold (Fig. 0). These extra states do not 



E 

MG spinon continuum MG spinon continuum 



Micro-spinon band 



VB excitations VB excitations 

I I I I 

QDM ground state QDM ground state 

FIG. 6: A cartoon of the spectra for the limiting cases where 
the Klein ratio is large (left) and small (right). In the small 
Klein ratio case (right), the lowest excited states are described 
by extended spinon wavefunctions. These scattering states 
are present even in the large ratio case (left) but here the 
lowest excitations, which we call "microscopic" spinons, are 
described by a wavefunction having peaks at the reference 
sites and decaying on the chains. The decay rate can be made 
arbitrarily fast by tuning the Klein ratio. The bandwidth of 
these localized states depends on the decoration and is zero 
in the infinite decoration limit. Details are given in Appendix 

m 



arise when we consider the decorated square lattice. Con- 
sider perturbing a decorated square lattice Klein model 
with: 

5H = J s"i ■ sj + 

m 

v^2((si ■ s bl )(s 3 ■ s b3 ) + (si • s ai )(s 3 ■ s 03 )J 

□ 

(3.19) 

where the first term is a nearest neighbor interaction 
and the second sum is over all elementary plaquettes, 
the spins 1234 labelling the 4 sites of a square plaquette 
in clockwise order and, as before (see Fig. [SJ), the labels 
ai and bi denoting the first counterclockwise and first 
clockwise neighbor of spin i. By arguments similar to 
the honeycomb case, we obtain an effective Hamiltonian: 

H af 3 = -Jx i{N+1) U a p + vn fl , a 5 a/ 3 

= -tu a0 + vn f i ia 6 a p + 0{vx i(N+1) + tx 2N ) 

(3.20) 

Here Uij is a matrix that is 1 if states \i) and |j) differ by 
the (minimal) length 4 loop and zero otherwise. There- 
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fore, we realize the physics of the square lattice QDM. 
Note that without the decoration, the error would be or- 
der x 2 = (1/2), as opposed to the bare honeycomb case 
where the error is order x 4 = (1/4). 




FIG. 7: The Klein model for the square lattice, with periodic 
boundary conditions, permits ground states which are not 
in manifold spanned by nearest-neighbor valence bond states 
such as this one. Here the thin lines form the lattice and 
the thick lines denote singlet pairings. Note that the Klein 
condition is satisfied at every lattice site. 




FIG. 8: The Klein model for the triangular lattice admits 
many nontrivial non-dimer ground states, such as this one. 
The thin lines show the lattice and the thick lines denote 
singlet pairings. The dots represent free spins. Note that the 
Klein condition is satisfied at every lattice site. 



IV. OTHER VALENCE BOND PHASES IN d = 2 

AND d = 3 

We now sketch the application of our strategy to obtain 
the rest of dimer model physics, including RVB phases. 
The points to be made concern the choice of lattices and 
phase conventions. 



A. Non-bipartite lattices in d = 2 

The simplest d = 2 non-bipartite lattice is the trian- 
gular lattice. Numerical evidence suggests that (nearest- 
neighbor) valence bond states are linearly independent. 25 
As with the square and honeycomb lattices, we may ori- 
ent bonds so that states differing by a (minimal) length 
4 loop, have positive overlap^ As with the square lat- 
tice, the Klein model admits non-valence bond ground 
states, though the problem is more serious with the tri- 
angular lattice (see Fig. [SJ. Decoration eliminates these 
possibilities by removing the triangular nearest-neighbor 
structures. Applying our strategy to the decorated tri- 
angular lattice allows us to reproduce the physics of the 
triangular lattice QDM, including its RVB phaseAi. By 
calculations similar to Appendix [5] one may show that 
spinon excitations are gapped. For the triangular lattice, 
it is known that collective excitations within the valence 
bond manifold are also gappecUi^, a conclusion which 
will remain valid for the decorated case. Therefore, we 
have constructed a model that shows a stable, SU(2)- 
invariant RVB liquid phase. 

Another non-bipartite lattice is the pentagonal 
lattice^ shown in Fig. [§] There are currently no for- 
mal proofs for the pentagonal lattice regarding the is- 
sues of linear independence of (nearest-neighbor) valence 
bond states and whether the set of these states spans 
the ground state space of the corresponding Klein model. 
However, it was explained in Ref.T6 that the most impor- 
tant ingredients of their proofs for the honeycomb lattice 
are its relatively low coordination number (3); relatively 
large minimum loop size (6); and the absence of triangu- 
lar structures in the lattice. The pentagonal lattice, has 
sites of coordination 3 and 4, minimum loop size 8, and 
no triangular structures, suggesting that the arguments 
may be adapted to this lattice. As before, it is possible 
to orient bonds so that the overlap of states differing by a 
(minimal) length 8 loop always has the same sign. While 
in the square, triangle, and honeycomb cases, the sign 
of the overlap is a matter of convention (which we chose 
as positive), for the pentagonal lattice, only the negative 
sign is possible. In fact, in Appendix 1X1 it is shown that 
using the fermionic convention, one may always obtain 
the negative sign independent of lattice details. There- 
fore, to generate the QDM kinetic energy, we must per- 
turb the Klein model with a ferromagnetic interaction. 
From a perturbation of the form: 

5H = — J Si ■ Sj + 

(ii) 

«y^((si • s bl )(s 3 ■ s ba )(s 5 ■ Sb 5 )(sV • s b7 ) 

o 

+(si • s 0l )(s 3 • s aa )(s 5 ■ s a5 )(s 7 ■ s ar )j (4.1) 

where the first term is over all nearest-neighbor spins 
and the second term is over elementary plaquettes {fli 
and bi, once again, denoting the first counterclockwise 
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and first clockwise neighbors of spin i), we may obtain 
the quantum dimer Hamiltonian for the decorated pen- 
tagonal lattice: 

H aP = -tO af} + vn flta 5 aP + 0(vx 8 ( N+ V + tx 2N ) 

(4.2) 

Here Oy is a matrix that is 1 if states \i) and \j) dif- 
fer by the (minimal) length 8 loop and zero otherwise. 
Therefore, we realize the physics of the pentagonal lat- 
tice QDM. We have checked that in the classical limit, the 
dimer-dimer correlations decay exponentially and in Ap- 
pendix El we present numerical evidence that monomers 
are deconfined. As both features also hold at the RK 
point, we may repeat the arguments described in Ref. ^] 
for the triangular lattice to conclude that the pentago- 
nal lattice QDM also shows an RVB liquid phase, a result 
which may be transcribed into spin language as discussed 
above. 




FIG. 9: The pentagonal (Sutherland-Shastry) lattice. It is 
the dual of the Shastry-Sutherland lattice, which is indicated 
by the dashed lines. 



B. Non-bipartite lattices in d = 3 

The face-centered cubic (FCC) is a three-dimensional 
non-bipartite Bravais lattice with each site having 12 
nearest- neighbors. The undecorated lattice has triangu- 
lar structures involving two neighboring facial sites and 
the two corners which are their common neighbors. This 
will lead to non-valence bond ground states in the FCC 
Klein model. Decoration eliminates the triangular struc- 
ture and hence this type of pathology. The shortest reso- 
nance loops are length 4. In the fermion sign convention, 
these loops come with negative sign. A perturbation con- 
sisting of a ferromagnetic exchange and 4-spin resonance 
interaction will reproduce the QDM results for the dec- 
orated lattice. We expect the resulting model to show a 
Zi RVB phase near its RK poinliiS. 



C. Bipartite lattices in d — 3 

For the diamond lattice we pursue the same strategy as 
above. The properties of the diamond lattice we require 
are the following. 

The diamond lattice is bipartite, has coordination four, 
and the shortest resonance loops are of length six. The 
Klein Hamiltonian again has nearest-neighbour dimer 
coverings as ground states, although, as for the case of 
the pentagonal lattice, no theorem exists excluding other 
ground states. It is likely that extra states may be ex- 
cluded by decorating the lattice. The number of dimer 
ground states, n gs , is exponentially large in the number 
of sites, N, but it is not known exactly. Defining the 
ground state entropy per site as S = (1/N)hm gs , an 
accurate series expansion by Nagle^ yields S « 0.265. 

We now try to mimic an RK quantum dimer model 
for the diamond lattice. As for the case of the honey- 
comb lattice, we do this by adding a nearest-neighbour 
exchange term to induce a kinetic term and in addition, 
a ring term to generate a potential term. We then ex- 
pect the resulting model to exhibit, near the effective 
RK point, a U(l) RVB liquid phase with algebraically 
decaying correlations as well as gapless photonic gauge 
excitations, as discussed in detail in Ref. [01 

This liquid phase will give way, upon making the po- 
tential term more attractive, to a columnar-type solid. 
For an increasingly repulsive potential, the scenario of 
Cantor deconfinement predicted for the two-dimensional 
case is simplified. Technically, there are no relevant lock- 
in terms in three dimensions so that the deconfined region 
simply acquires an increasing amount of U(l) flux as v/t 
is increased through the RK point; finally, a staggered 
solid, with the maximal amount of U(l) flux allowed by 
microscopic constraints ji2iiiLi£ is reached. We cannot say 
whether this will happen continuously or via a first order 
transition. 



V. DYNAMICAL SELECTION OF GAUGE 
STRUCTURES: PYROCHLORE LATTICE 

We construct a Klein model with Zi order and a 
Kivelson-Klein model with U(l) order which takes advan- 
tage of the bipartiteness of the dual lattice. This nicely il- 
lustrates the dynamical selection of the low-energy gauge 
structure present in topological phases. 

A. The Klein model 

The undecorated pyrochlore lattice (Fig. I1L)[1 does not 
lend itself straightforwardly as a starting point for dimer 
models obtained via the Klein route because its basic 
building block, the tetrahedron, supports more dimer 
coverings than linearly independent singlet states. By 
sufficiently decorating the lattice, the orthogonality prob- 
lem is solved. The shortest resonance loops are length 6 
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FIG. 10: The pyrochlore lattice, a network of corner-sharing 
tetrahedra 

and the fermionic convention may be used to make the 
minimal overlaps come with negative sign. Perturbing 
the decorated pyrochlore Klein model with a ferromag- 
netic nearest-neighbor interaction and 6-spin ring inter- 
action, we obtain an effective Hamiltonian mimicking the 
pyrochlore lattice QDM, which includes a RVB phase. 

B. The Kivelson-Klein model 

A modified version of the Klein Hamiltonian, which we 
will refer to as the Kivelson-Klein Hamiltonian, may be 
used to produce a model displaying a (7(1) RVB phased 
Its Hamiltonian is of the same form as Eq. 13.11 but the 
definition of M(i) is changed. The projection now acts 
not on a site and its nearest neighbours but instead on 
the four sites of a tetrahedron: 

tet 

We note that the simple decoration trick described above 
cannot as usefully be applied to the Kivelson-Klein 
model, as here increasing the number of sites in the tetra- 
hedron does not lead to an increase in the number of 
dimers required. 

1. Ground states of the Kivelson-Klein model 

Evidently, each state in which each tetrahedron con- 
tains at least one singlet bond is a ground state of the 
above Hamiltonian, Eq. 15.11 How can this be related to 
dimer coverings of the pyrochlore lattice? 

First, note that (i) the number of tetrahedra equals 
twice the number of sites, N t — 2N, and that (ii) the 
number of hardcore dimers, Nd < N/2, as a dimer in- 
volves two sites. From this it follows that Nd < N t , 
the equality sign holding for hardcore dimer coverings. 
However, the requirement of having at least one dimer 
in each tetrahedron gives (iii) Nd > Nt. We therefore 
see that (ii) and (iii) imply that those hardcore dimer 
coverings of the pyrochlore lattice for which each tetra- 
hedron contains exactly one dimer are ground states of 
the Kivelson-Klein Hamiltonian. 

The ensemble of these states maps onto the ground 
states of the six-vertex model on the diamond lattice, 



or equivalently onto the ground states of the pyrochlore 
Ising antiferromagnet. This can be seen as follows. First, 
note that the lattice of tetrahedra defined by the py- 
rochlore lattice is the bipartite diamond lattice. One 
diamond sublattice sits at the centre of the 'up' tetra- 
hedra, the other one at the centres of the 'down' tetrahe- 
dra. Now let us define the Ising spins on the pyrochlore 
lattice as follows. For an up (down) tetrahedra, the pair 
of spins at the two ends of a dimer point up (down) , and 
the other pair points down (up). This defines a one-to- 
one mapping of dimer to Ising states; crucially, on each 
tetrahedron, two spins point up and two point down, thus 
putting the tetrahedron into an Ising ground state. (The 
mapping to the six-vertex model on the diamond lattice 
proceeds by calling an up (down) spin an arrow pointing 
from the centre of an up (down) tetrahedron to a down 
(up) tetrahedron). 



2. Ground state correlations 

The total number of ground states (assuming that 
there are none in addition to the abovementioned dimer 
states) gives rise to an extensive ground state entropy 
well-approximated by the Pauling entropy Sp au u ng — 
(l/2)ln(3/2). 

Using the mapping to an Ising magnet, it is straight- 
forward to calculate the correlator between singlet bonds 
averaged over the ground state manifold. To do this, note 
that each of the six dimer positions on a bond of a given 
tetrahedron corresponds to an Ising ground state of that 
tetrahedron. In turn, the corresponding vertex of the six- 
vertex model describes a net flux, the direction of which 
is given as follows. Consider a cube circumscribing the 
tetrahedron in question, so that the bonds of the tetrahe- 
dron are face diagonals of the circumscribing cube. The 
direction of the flux (i.e., the average direction of the four 
arrows of the given vertex) now points from the centre of 
the cube through the midpoint of the face of which the 
bond occupied by the dimer resides. 

Using the theory developed in Ref. one can read off 
that the dimer correlations are simply dipolar. Briefly, 
this follows from the observation that, upon coarse grain- 
ing, the smaller the coarse-grained flux, the more mi- 
crostates (prior to coarse-graining) correspond it. Mod- 
elling this by an effective quadratic weight on the flux 
configurations leads to simple magnetostatics. 

For example, the connected correlator between a pair 
of dimers located on the top of an up tetrahedra, sepa- 
rated by a vector r which makes angle 9 with the z axis, 
is proportional to the dipolar form (3cos6> 2 — l)/r 3 . 

Finally, this model can again in principle be "Rokhsar- 
Kivelsonized" , i.e. by adding appropriate perturbations 
to Eq. 15.11 we may generate an effective Hamiltonian 
which acts on the Kivelson-Klein ground state manifold 
in a manner similar to Eq. 12.11 on the space of dimer 
coverings. We do not do this here for a number of rea- 
sons. Firstly, the expected phase diagram has the same 
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topology as that discussed for the diamond lattice in the 
previous section, so no new phases are obtained. Sec- 
ondly, the dimer dynamics is rather messy. The short- 
est resonance loop now involves six dimers straddling a 
hexagonal loop of the pyrochlore lattice - and there are 
several symmetry-inequivalent loops of this type. (This 
is reminiscent of the - exactly soluble - kagome dimer 
model proposed in Ref.l37l) In addition, for a simple per- 
turbing nearest-neighbour exchange, the resonance term 
vanishes in the leading order of the overlap expansion, so 
that a more complex perturbation is needed. 



VI. DISCUSSION AND OUTLOOK 

We have presented spin-1/2 Heisenberg Hamiltonians 
that realize a large class of valence bond phases. In par- 
ticular they realize Zi RVB phases in d = 2 and d = 3, 
the U(l) RVB phase in d = 3 and the Cantor decon- 
fined region in d = 2. These phases have previously been 
shown to exist in quantum dimer models with dimers 
standing in for valence bonds. In this paper we have 
constructed Klein models that exhibit ground state man- 
ifolds spanned by nearest neighbor valence bond states 
and then perturbed them to to realize quantum dimer 
models within these manifolds. This perturbation is done 
within the framework provided by the overlap expansion, 
made arbitrarily accurate by a decoration procedure that 
we have introduced. 

The decoration has the effect of expanding the length 
scale on which the Hamiltonian acts directly. However in 
order to stabilize the phases in quantum dimer models, 
we do not need to go to infinite decoration - it would be 
enough to suppress subleading terms sufficiently. In this 
fashion we obtain spin models with interactions of finite 
ranged While large decorations are needed to realize the 
simplest dimer models under analytic control, there is lit- 
tle reason to doubt that the extent of decoration could 
be reduced drastically, and even eliminated altogether 
on some of the lattices without sacrificing the various 
phases of interest. The subleading interactions will not 
neccessarily uniformly tend to subvert such phases and 
it is also possible to add other terms that would stabilize 
them. Showing how this can work is an obvious task for 
the future, as is the construction of mathematically rig- 
orous proofs for various statements in this paper that are 
made by appealing to a small parameter. One promising 
approach is that of Ref. |H, which appeared at the same 
time as our work, where the phases of the d — 2 bipartite 
quantum dimer model are realized by perturbing a Klein 
model with ring exchange terms. 

We emphasize that our central result has been the 
demonstration that spin liquid phases can be realized 
in SU(2)-invariant models. However, the actual mod- 
els we have given involve rather complicated geometries 
and Hamiltonians without direct experimental relevance. 
While it may be possible to engineer highly decorated 
lattices, a more important task perhaps, now that the 



question of principle is settled, is to refocus on studying 
much simpler Hamiltonians. Our insistence on a specified 
form of the wavefunction (containing nearest neighbor va- 
lence bonds alone) has led to fairly complicated Hamil- 
tonians but simpler Hamiltonians can exhibit the same 
phases with more elaborate ground state wavefunctions. 
Indeed, the situation for the simplest lattice under con- 
sideration in this paper, the (undecorated) honeycomb 
lattice, looks quite promising. Previous exact diagonali- 
sations of a J\ — J2 — J3 Heisenberg model on this lattice^ 
have clearly demonstrated the existence of the staggered 
VBS. The current data appears not inconsistent with a 
scenario in which the magnet leaves the staggered phase 
but never reaches the fluxless plaquette phase. We are 
optimistic that it will be possible to realize the physics 
discussed in this paper in such simpler models. 



Acknowledgements 

We would like to thank Chris Henley for pointing out 
to us the pentagonal lattice discussed above. We are very 
grateful to Steve Kivelson for formulating the generali- 
sation of the Klein model on the pyrochlore lattice. We 
would like to thank Matt Hastings for a highly stimulat- 
ing discussion which initiated our thoughts regarding the 
decoration procedure. We further thank Lincoln Chayes, 
David Huse and Achim Rosch for useful discussions. RM 
is grateful to the Aspen Center for Physics, where parts 
of this work were undertaken. This work was in part 
supported by the Ministere de la Recherche et des Nou- 
velles Technologies with an ACI grant. SLS would like 
to acknowledge support by the NSF (DMR-9978074 and 
0213706) and the David and Lucile Packard Foundation. 



APPENDIX A: SIGN CONVENTIONS IN THE 
OVERLAP MATRIX 

In the Rokhsar-Kivelson derivation of the dimer model, 
the overlap matrix between the different dimer cover- 
ings plays a crucial role. Whether or not the RK point 
corresponds to an equal-amplitude superposition of all 
dimer coverings depends on the question of whether the 
dimer Hamiltonian can be turned into a form where all 
off-diagonal matrix elements are negative. 

The leading effect of a perturbing nearest-neighbour 
exchange finally gives rise to a constant (for each type of 
loop) multiplying the overlap matrix (restricted to that 
type of loop). If the overlap matrix can be written as 
a term proportional to matrix with entries only or 1, 
the ground state is indeed given by an equal-amplitude 
superposition at the RK point. 

In the following, we demonstrate that a fcrmionic sign 
convention may be used to generate the negative sign, 
independent of lattice details, for models involving va- 
lence bond coverings of the lattice. We show that for 
the honeycomb and diamond lattices, we may obtain the 
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positive sign as well. We also present a convention for 
the Kivclson-Klcin model which gives the negative sign. 

1. Overlaps in the fermionic convention 

A general convention for overlaps can be obtained by 
employing the so-called fermionic convention, where a 
valence bond between sites a and b is generated via op- 
erators, such as the one placing a fermion with spin up 
on site a: c^. The singlet bond is then defined as: 

IN) = <4lo) - ^iAi + 4 T 4j|o) • (Ai) 

Here |0) is the vacuum state with no fermions present. 
Note that 

dl = 4 a , 

and that these operators, being bilinear in fermions, com- 
mute unless they have exactly one site in common. This 
means that for constructing a valence-bond covering, the 
order in which the bonds are generated is inconsequen- 
tial. 

A loop in the transition graph involving sites h, i, j, 
and k will lead to the following type of expression in the 
overlap matrix element calculation: 

(o\d ab ---d ij 44- • -4 6 |o) 

= ~<0|da6-"4*-"46l0> ( A2 ) 

By induction, a loop in the transition graph involv- 
ing L dimers in each configuration leads to a factor of 
(— l/2) i_1 , independently of any further details of the 
lattice. 



2. Honeycomb and diamond lattices 

The following approach works for both honeycomb and 
diamond lattice, both of which have a shortest resonance 
loop of length six, as in the original benzene picture. 

As these lattices are bipartite, we can orient each bond 
to point from one sublattice (A) to the other (B), so that 
a singlet between sites a and b of sublattices A and B, 
respectively, has the wavefunction 

\(ab)) = [\UU)-\ laU)]/V2 . (A3) 

The two singlet coverings of the benzene loop now have 
wavefunctions 

|1> ee \(ab)(cd)(ef)) ; \2) = \(bc)(de)(fa)) 

from which it follows that 

<1|2> = +1/4 (A4) 



for any hexagonal plaquette. 

It is in fact also possible to choose 

(1|2> = -1/4 (A5) 

for the honeycomb lattice. This can, for example, be 
achieved by choosing any fixed hardcore dimer covering 
of the triangular lattice which is dual to the honeycomb 
lattice under consideration. One then multiplies each 
valence bond state of the honeycomb lattice by (— 1)™ X , 
where n x is the number of valence bonds which cross 
dimers of the triangular dimer covering. This generates 
the desired effect. 



3. Other bipartite lattices 

The above construction for generating uniform overlap 
matrix elements can be generalised to any bipartite lat- 
tice. By orienting the bonds from one sublattice to the 
other, one always obtains an overlap which is positive; 
its size is the simple product over the individual loops 
involved in the transition graph. 

(k\l) > ; (A6) 

indeed, the positive overlap holds true for any value of 
the loop length, and therefore for an arbitrary pair of 
valence bond coverings |fc), \l). 

4. Kivelson-Klein model on pyrochlore lattice 

Here we first need to establish the possible resonance 
loops. These involve six dimers on a cluster of twelve 
sites arranged as follows. 

Six of the sites sit on a hexagonal ring on the py- 
rochlore lattice; each link of this hexagonal ring belongs 
to a different tetrahedron. Each of these six tetrahcdra 
contains one dimer linking a site on the hexagonal ring 
with a site off the hexagonal ring. As there is a choice 
of 2 such off-sites per tetrahedron, the total number of 
shortest resonance loops corresponding to each hexagonal 
ring equals 2 6 = 64. Not all of these loops are symmetry 
equivalent. 

These loops all involve moving six dimers. Hence, their 
overlap in the Fermion convention is given by —1/32. 

APPENDIX B: SPINON GAP FOR THE 
DECORATED HONEYCOMB LATTICE 

A stable RVB liquid phase requires not only certain 
properties of the ground state wavefunction but also that 
the excitation spectrum has a lower bound. In this sec- 
tion, we argue that the nearest neighbor valence bond 
ground states are separated by a finite gap from the ex- 
cited states for the case of the decorated honeycomb lat- 
tice Klein model. 
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To this end observe that in the highly decorated limit 
we are connecting a set of Majumdar-Ghosh 7 (MG) 
chains into a two-dimensional network. The excitations 
of the chains themselves are well studied: these are 
spinon defects between the two different dimerizations 
and there is considerable analytical^^ and numerical^! 
evidence that they are gapped. In putting the chains 
together we need to ask if the intersections lead to the 
emergence of states below the one-spinon continuum on 
the chains that can fill in the gap. In the infinite decora- 
tion limit, it is sufficient to consider a single intersection: 
thus we look for bound states localized near a site of the 
original honeycomb lattice where three MG chains would 
cross (see Fig. ITT|) . 




FIG. 11: This figure depicts a spinon at the crossing of the 
edges of a decorated honeycomb lattice. 



We consider a variational wavefunction describing a 
single spin at the crossing of three semi-infinite MG 
chains. 



>0|00) e 



(Bl) 



where |+)o denotes an up spin at location and |00) ee 
denotes that "everything else" is dimerized into singlet 
(00) pairings. We consider the action of the Klein Hamil- 
tonian l|3.1|l on this state. Observe that JV(0)> which 
involves only site and its three neighbors, is the only 
projector that does not destroy our trial function. A 
standard calculation gives: 
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iV(0)|V»>0 = C cr[-J Mo+4(JV>2a- 
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acting on a spinon state produces the original state along 
with states where the spin has hopped two sites. We also 
obtain terms which are orthogonal to all spinon states, 
designated here by the expressions \a) and \ j3). This mo- 
tivates the bound state trial function: 
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where y is a variational parameter less than unity. In 
calculating the expectation of the energy, we need to con- 
tend with the non-orthogonality of the spinon states: 



(B5) 



for sites i and j on the same (different) chain(s). A te- 
dious but standard calculation^ gives the expectation 
value for the energy of the trial state l|B4|l : 
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where \ip)i — |+)i|00) ee describes an up spin at location 
i with everything else paired into singlets and C cr sets 
the energy scale for violating the Klein condition at the 
cross. We can also compute the effect of operator (|3.1f) 
on the state \ip)i, where i is a location along a chain; in 
this case, only contributes: 



= a 
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where a — C c h/C cr . 

Fig. El sriows a graph of this expression for several val- 
ues of a. For small values of a, the spinon wavefunction 
is an extended scattering state while for large values, the 
wavefunction is localized at the cross. 

Our analysis has been for the infinite decoration limit. 
For large, but finite, decoration, the extended spinon 
states obtained for small a may be interpreted as MG 
spinons scattering at the vertices. In this limit, the nat- 
ural extension of the localized states obtained for large 
a will involve the wavefunction having peaks at the ref- 
erence sites and decaying on the chains. There will be a 
band of such localized states below the scattering states. 
We may estimate the bandwidth by considering a varia- 
tional wavefunction where the spinon resides only on the 
reference sites: 



Here C c h is the energy scale for violating the Klein con- 
dition on the chain. We see that the Klein Hamiltonian 
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FIG. 12: This plot gives the energy expectation as a func- 
tion of the variational parameter y for different values of 
a = C c h/C C r, which is ratio of Klein scales; the energy is 
measured in units of C cr - For large values of this ratio, the 
minimum occurs for y — which corresponds to spinons local- 
ized at the crossing. For small values, the minimum occurs at 
y = 1, which is an extended spinon state. At a value slightly 
greater than a = 1, there is a first order "phase transition" 
between these limits. The important feature is that for any 
nonzero a, there is an energy gap between spinon states and 
the valence bond states, which have zero energy. 



where \ri) denotes a wavefunction of the form (|B1(I for 
the spinon at lattice site n. The variational energy of 
this trial state, to leading order in a large N expansion, 
where N is the number of sites inserted between reference 
sites in the decoration, is given by: 



E — C c 
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where x = -h=. We see that the band becomes more 
narrow as the decoration is increased. 

In our analysis so far, we have considered states where 
the defects are always an even number of sites away from 
the reference sites. There is another family of spinon 
states corresponding to the defects being located on the 
odd sites. Referring to Fig. 1131 the Klein operator per- 
mits the spin at la to hop only to the site lb, which 
is connected by a dimer to the origin. Therefore, in the 
large decoration limit, this is equivalent to the MG chain, 
which we know is gapped. 



APPENDIX C: RESULTS FOR THE 
PENTAGONAL LATTICE 

At infinite temperatures, where "infinite" means a 
temperature that is small compared with the excita- 
tion gap of the Klein Hamiltonian but much larger than 




FIG. 13: This shows a representative of the family of states 
where a spinon is an odd distance from the origin. In this 
configuration, the Klein operator may hop the spinon only 
onto the 6-chain. The configuration where the origin forms a 
singlet with lc is essentially orthogonal to the given configu- 
ration for large decoration. 
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FIG. 14: The monomer-monomer correlation function for a 
100x100 pentagonal lattice with periodic boundary condi- 
tions. The two curves are cuts along the x and x + y di- 
rections. The distance r refers to the distance between unit 
cells. In computing these correlation functions, we take the 
two test monomers to be on the same sublattice (the pentag- 
onal lattice is a cubic lattice with a six point basis). Each 
data point is an average over N = 10 6 data points and the 
noise seen in the plot is of the order of Monte Carlo noise 

i/Vn ~ io" 3 ~ 0.1%. 



the energy scales of the quantum dimer Hamiltonian, 
the dimers are described classically, i.e. thermodynamic 
quantities are computed as equal-weight averages over all 
dimer states. The number of dimer states grows expo- 
nentially with lattice size. This number may be com- 
puted using the method of Kasteleyn4iii£ The results 
are shown in tabled The method also yields the entropy 
per site in the thermodynamic limit. 



5 = 0.168608.. 



(CI) 



The striking feature of table [i] is that even fairly small 
systems have an enormous number of dimer coverings 
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TABLE I: Table of the number of dimer configurations for an 
L x xL y size pentagonal lattice. Here L x and L y refer to the 
underlying square lattice; the pentagonal lattice is a square 
lattice with a 6 point basis. We have assumed periodic bound- 
ary conditions in this calculation. 



so numerical studies of large systems require Monte- 
Carlo simulations. Fig. ^] is a Monte-Carlo calcula- 
tion of monomer-monomer correlation functions C (r) for 
the pentagonal lattice using the algorithm of Sandvik 43 . 
C(r) is defined as the number nd(r) of dimcr coverings 
given a pair of test monomers separated by distance r 
divided by n^(l). The simulation shows monomer de- 
confinement at infinite temperatures, as opposed to the 
square and honeycomb lattices which show logarithmic 
confinement. This indicates a liquid phase at high tem- 
peratures. The RK point also has deconfined monomers 
and the same dimer correlations as the high tempera- 
ture phase, which strongly suggests that the RK point is 
part of a zero temperature liquid phase which connects 
continuously to the high temperature liquid. 
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